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ABSTRACT 

We present a detailed kinematical analysis of the young compact hourglass-shaped planetary 
nebula Hb 12. We performed optical imaging and longslit spectroscopy of Hb 12 using the 
Manchester echelle spectrometer with the 2.1m San Pedro Martir telescope. We reveal, for the 
first time, the presence of end caps (or knots) aligned with the bipolar lobes of the planetary 
nebula shell in a deep [Nil] A6584 image of Hb 12. We measured from our spectroscopy 
radial velocities of ~ 120 kms^' for these knots. 

We have derived the inclination angle of the hourglass shaped nebular shell to be ~ 65° 
to the l ine of sight. It h as been suggested that Hb 12's central star system is an eclipsing 
binary (iHsia et al.ll2006l) which would imply a binary inclination of at least 80°. However, if 
the central binary has been the major shaping influence on the nebula then both nebula and 
binary would be expected to share a common inclination angle. 

Finally, we report the discovery of high-velocity knots with Hubble-type velocities, close 
to the core of Hb 12, observed in Ha and oriented in the same direction as the end caps. 
Very different velocities and kinematical ages were calculated for the outer and inner knots 
showing that they may originate from different outburst events. 

Key words: circumstellar matter - stars: kinematics - stars: mass-loss - ISM: jets and out- 
flows - planetary nebulae: individual: Hb 12. 



1 INTRODUCTION 

Understanding the wide variety of planetary nebulae (PNe) mor- 
phologies has puzzled the stellar evolution community for many 
years. The Generalized Interacting Stellar Winds (GISW) model 
has long been t he widely accepted model for the production of non - 
spherical PNe i lKahn & Westlll985l:lBalick, Preston & Ickdll987h . 
It involves a strong equatorial density enhancement in the red giant 
wind around the progenitor, into which runs the fast wind from the 
central white dwarf, which causes the fast wind to progress faster 
in the polar direction, thus giving the nebula its bipo lar shape. 

However, in recent years, Sahai & Trauger ( 1998) and 
ISoker & Rappapor^ 42000) have argued that extremely elongated 
nebulae and bipolar shells with very narrow waists could only be 
shaped by high-speed collimated outflows or jets rather than pre- 
existing equatorial density enhancements. 

* E-mail: neil.vaytet@cea.fr 



With the arrival of the Hubble Space Telescope (HST) and 
state of the art CCD instruments in the mid 1990s, high-resolution 
images of PNe also revealed the widespread presence of mi- 
crostruc tures such as bi- or multi-polar jets, knots, bubbles and fil- 
aments ( [Harrington & Borkows 020061) . 

Two major lines of investigation emerged into the origin of 
such microstructures. Single-star models attributed the collima- 
tion of outflows to strong magnetic fields which are fuelled by 
differential rotation in the stellar envelope. The claim was sup- 
ported by three-dimensional magneto-hydrodynamic simulations 
which were successful in creating an impressive range of PN mor- 
phologies including multi-lobes and point s ymmetric PNe, which 
could not be accounted for by GISW theory (Dwarkadas & Balicy 
1998; Garci'a-Segura et al. 1999; Garci'a-Segura & Lopez 200d). 
Recently, toroidal magnetic fiel ds have been directly obs erved in 
the proto-P N W43A ( Vlemmings, Diamond & Imai"2006l) and the 
PN K3-35 jMiranda e t al. 2001; Gomez et al. 2009) through maser 
measurements, and in PNe NGC 7027, 6537, 6302 and proto- 
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PN CRL 2688 th rough polarisation measurements ( lGreavesll2003 : 
ISabin et alj2007b . proving the existence of such fields. However, it 
has been argued that magneto-rotational launch, the mechanism for 
launching collimated jets with powerful m agnetic fields, requires 
more energy than a sing le star can provide ( Nordhaus et alJllOOTl : 
lFrankl2006| - fSoke JI2004I) . The star would initially be able to launch 
a jet but would not be able to sustain such magnetic fields for long 
periods of time, making the model inconsistent with PN evolution. 

Collimated outflows are usually associated with quasars or 
X-ray binaries, where the accretion disc plays an essential role 
in feeding mass, en ergy and angular momentum into the jet. 
ISoker & Liviol ( 1 19941) hence claimed that a stellar accretion disc 
must also be required to produce such outflows in PNe, thus infer- 
ring the presence of a binary companion. The outflow from a star 
in a close-binary system will also inevitably undergo a common- 
envelope (CE) stage during which the secondary is engulfed by 
the primary's expanding atmosphere into which it deposits energy 
and angul ar momentum through friction, g iving the outflow a bipo- 
lar sh ape jLlovd, O'Brien & Bodelll997l : IPorter. O'Brien & Bodd 
Il998h . 

In very recent years, binary systems have been considered to 
be essential to the formation of very elongated PNe. An ac cretion 
disc is a very efficient source of collimation dSokeJbOOlh and a 
companion star can spin-up the primary to help it sustain power- 
ful magnetic fields via angular momentum convervation. Obser- 
vational evidence of precessing jets in some PNe al so implies the 
prese nce of binary central systems jLopez. Vazquez & RodriguezI 
19951). Finally, recent st ellar populations statistics carried out by 



Moe & De Marcol ( |2006|) showed that if both single stars and bina- 



ries evolve to PNe, the number of observed PNe is much lower than 
expected. The numbers are much closer if only binary stars are al- 
lowed to evolve to PNe, suggesting that binary systems are in fact 
essential to forming all PNe. However, currently only a very small 
fraction of PNe have confirmed close-binary central stars, likely to 
have produced accretion discs (De Marco 2009). Esti mations sug - 
gest that 10%-15% of all PNe harbour close binaries (lBondll2000li . 
and another ~ 10% PNe with wider binaries ([CiarduUo et"al .119991 : 
see also the full discussion bv lDe Marcdl2009l) ! 

The CE phase and jet collimation by accretion discs predict 
that a bipolar PN should be extended in the direction normal to 
the binary orbital plane, and also imply that any PN with a cen- 
tral close-binary should be non-spherical. Little observational work 
has been carried out on PNe harbouring close-binary central stars. 
We started in 2004 an observational campaign in order to obtain 
morpho-kinematical models of all PNe with confirmed close-binary 
central stars. We uncovered for the first time a direct link between 
PN shell and binary syst em after a common i nclination angle was 
found in Abell 63 (A63) I Mitchell et a"i]|2007h . 

We present in this paper a kinematical study of the young PN 
Hubble 12 (Hb 12, R.A. 23:26:15, Dec. +58:10:55 [J2000.0]); an 
extremely well studied compact nebula 10" in diameter) which 
posse sses a bipolar hourglass shell morphology ( Sahai & TraugeJ 



19981) . It lies at a distance of 2 .24 kpc Icahn, Kaler & St angheUinil 



1992h . iMiranda & Solj fl989l) performed a kinematic analysis of 
the nebula and found the inner shell to be expanding at ~ 16km s^ 
and faint extended bipolar lobes. They also detected the presence 
of a faint, extended halo in Ha . 

,Hsia. Ip & Li ( 2006) detected periodic photometric variabil- 
ity in the optical lightcurve of the central star and attributed 
it to an eclipsing binary system. However, doubt was cast on 
this interpretation of the o bserve d photometric variability by 
iDe Marco. Hillwig & SmithI | |2008|) . suggesting that there might 



not even be an eclipsing binary or that at least the binary param- 
eters suggested were wrong. 

Our aim is to construct a three-dimensional kinematical model 
of Hb 12, through the use of longslit spectroscopy, in order to deter- 
mine its morphology and orientation. We will then be able to judge 
whether it is consistent with an eclipsing binary central star or not 
(assuming the binary is the main influence on the PN morphology). 
We also reveal for the first time direct observations of high- velocity 
knots in the system, both at a large distance and close to the central 
star, emitting in the optical; signs of strong outflow activity. 



2 OBSERVATIONS 
2.1 Imaging 

A deep 900 s narrowband [Nil] A6584 imag e of Hb 12 was taken 
using the Manchester echelle spectrometer jMeabum et al.ll2003h 
combined with the 2.1-m, f/8 San Pedro Martir telescope (MES- 
SPM) on 30* June 2000; this is shown in Fig.^ The TEK-1 CCD 
was used with 1x1 binning (= 0.31"/p ixel) and the seein g was 2". 
A high resolution HST image taken by ISahai & Trauged | |199S) IS 
also displayed in lower-right inset in Fig.[T] showing the hourglass 
shape of the core of Hb 12. 

The wide-field image unveils for the first time the presence 
of two knots on each side of the bright PN shell (indicated by the 
arrows) emitting in the low-ionisation [Nil] 16584 line, aligned 
with the polar axis of the PN. The lower left-inset shows the re- 
gion around the southern knot at a much higher contrast, in which 
three more much fainter knots (also indicated by arrows) begin 
to appear Such knots are often found around PNe and are usu- 
ally associated wi t h high velocity outflows (see 'Lope z et al.ll 19931 : 
iBalick et al.ll 19931 : lO'Connor et al. 2000; Mitchell et al.ll2()()7l: for 
example). In analogy to other hourglass PNe. lKwok & Hsial ( |2007|) 
predicted that outflow material should be present at a distance from 
the main nebula; our observations confirm this and are the first di- 
rect evidence for optical o utflow material around Hb 12. 

iKwok & Hsial ( |2007|) observe in their infrared HST image of 
Hb 12 equatorial features associated with a second outer hour- 
glass of a much bigger size and a wider opening angle, drawing a 
strong analogy to th e spectacular nested hourglasses PN Hen 2-104 
jCorradi et al.ll200ll) . They also compare it to MyCn 18 in which a 
nested hourglass was also found ( Sahai et al. 1999). The outer hour- 
glass in Hb 12, along with the 'eye' created by the intersection of 
inclined bipolar lobes (which strongly resembles th e central region 
of My Cn 18), are also vi sible in the H2 imag es of iHora & Lattej 
( Il996l) . However, as with lKwok & Hsial {2003), we do not see any 
signs of this second hourglass in our optical image. 



2.2 Longslit spectroscopy 

Longslit spectroscopy of Hb 12 was performed at high spatial 
and spectral resolution using the MES-SPM on the nights of 
10* November 1998 and 29* October 2000. Observations were 
taken through a narrow-band filter which isolates the Ha and 
[Nil] A6548, A6584 emission lines in the 87* echelle order For the 
1998 dataset, spectra were obtained using a TEK-1 CCD detector 
comprising 1024 x 1024 square pixels with sides measuring 24 jlm, 
using 2x2 binning (= 0.62"/pixel and 3.5 km s^Vpixel). A TH2K 
CCD was used for the 2000 dataset with 2048 x 2048 squai'e pix- 
els with sides measuring 14 ^um, using 3x3 binning (= 0.54" /pixel 
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Figure 1. A deep [Nil] A6584 image of Hb 12 taken with MES-SPM. North is up and east is left. The inset to the lower-right is an HSTAVFPC2 Ha enlargement 
of the core revealing the hourglass morphology of Hb 12 (HST aixhive / Saha i & Trauger 1998). The seven E-W slit positions are marked with dashed lines. 
The arrows indicate the positions of the knots. The inset to the lower-left shows the region around the southern knot, reveahng fainter knots, also indicated by 
arrows. 



and 3.1 km s /pixel). Integration times were 1200 s and the seeing 
was 2" for all spectral observations. 

A single slit with position angle PA = 86° was placed across 
the nebula at seven different positions (see Fig. [TJ; positions 2-4 
correspond to the data taken in 1998 whilst positions 1, 5, 6 and 7 
belong to the 2000 dataset. For both epochs, the slit was 150 /im 



wide (= 1.9" and 11 km s^') and 5' long. Slit positions 4 and 5 
intersect the south and north lobes, respectively, of the bright hour- 
glass core, and their [Nil] A6584 position-velocity (p-v) arrays are 
shown in Fig. (2] row 1, columns c and d, respectively. The Ha p-v 
arrays from slit positions 2 to 6 are shown in Fig. [3] Slit positions 1 
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and 7 go through the knots to the south (only the brighter one) and 
north of Hb 12, respectively. 

Data reduction was performed using the STARLINK software 
packages KAPPA, FIGARO and CCDPACK. All images and spec- 
tra were bias-corrected and cleaned of cosmic rays. The Ha + 
[Nil] spectra were wavelength calibrated against a ThAr emission 
line lamp. All velocities have been corrected to the heliocentric 
frame of reference. 



3 RESULTS 

3.1 Measuring the inclination of tlie hourglass 

There is no doubt about the core morphology of Hb 12, the hour- 
glass structure is clearly visible in the HST image inset in Fig.[T] 
However, the orientation of the shell is not clear. The common as- 
sumption for stellar outflows with binary central stars is that the 
orientation of the shell is close, if not identical, to t hat of the 
binary system (see iMitchell etaP 120071 for example). iHsia et aP 
( l2006l) suggested that the central star system of Hb 12 was a par- 
tially eclipsing binary after detecting periodic variability in their 
photometry, which would limit the binary's inclination angle to 
i'binai-y > 80°. Yet, on the HST images, the o rientation of the shel l 
appears to be closer to 50° — 70° , as shown bv lKwok & Hsi3 ^20071) 
who measured an inclination of /shell = 52° from fitting ellipses to 
co-axial rings observed in [Nil] HST images of Hb 12. Here we 
derive an inclination angle for the nebular shell based on the com- 
parison of model spectra derived from an hourglass shaped model 
nebula with our observed [Nil] spectra from slits 4 and 5 which 
cross the two lobes of the hourglass. 

We creat ed a synthetic hourglas s model of the shell using the 
SHAPE code jSteffen & LopezlboOfiO . The hourglass morphology 
of the PN shell was based on the HST imagery of Hb 12 (inset 
in Fig. [T] and the [Nil] images of iKwok & Hsial r2007h assuming 
axi-symmetry around the polar axis, and we gave it a Hubble-flow 
velocity law (proportional to dist ance from centre) com bined with a 
systemic velocity of —5 kms" 1 favung&Alle3ll996h (see Fig.[2j 
column a). In the HST image, the PN shell is visibly brighter on the 
western side; this asymmetry in brightness was reproduced in the 
model, with a 300% increase in brightness for the entire western 
half of the shell. The model was then convolved to the SPM image 
seeing of 2" (Fig. |2j column b), where the western brightening is 
now much more visible. We note here that extra bright emission is 
visible at the centre of the nebula in both the SPM and HST images, 
probably from the central star, which is not included in our model. 
However, this should not affect the following results since none of 

our slits positions go across this central region. 

As an aside, it is interesting to note here that lSahai & TraugeJ 

( Il998 ) found that the very core of Hb 12 (measuring ~ 0.6") it- 
self also shows signs of bipolarity with east-west inequalities in 
brightness, but unlike the PN shell, it is the eastern lobe of the core 
which is brighter than the western one; the origin of such anisotropy 
in bri ghtness is not kn own. Note that the structure observed by 
ISahar& TraugeJ ( Il998h in the core is large compared to the size 
of a close-binary system: the centres of the two lobes, assuming 
a distance of 2.24 kpc, are separated by ~ 350 AU. In addition to 
this, the elongated nature of the lobes indicates that they are most 
probably not observing a resolved stellar binary. 

Synthetic longslit p-v arrays were extracted from the model 
along slits matching positions 4 and 5 of the observations, con- 
volved to the velocity resolution of the MES and thermal broad- 
ening of the [Nil] A6584 line. Those were compared to the 



[Nil] A6584 p-v arrays obtained from the longslit spectroscopy 
(Fig. m row 1, columns c and d). They were compared to the 
[Nil] A6584 emission rather than Ha since Hb 12 is very bright 
in Ha and very little core structure is visible in the spectra at that 
wavelength, due to over-exposure. The Ha profiles were also con- 
taminated by scattering effects from the very bright core of Hb 12. 
In addition, the [Nil] A6584 line has an intrinsically narrower ther- 
mal width and tends to come from an outer 'skin' whereas the 
Ha emission tends to be distributed through the volume, which 
also broadens the profiles. 

Velocities of up to —40 kms^' and +30 kms^' are observed 
for positions 4 and 5, respectively, suggesting much higher ex- 
pansio n velocities than the 16 kms^' quoted by iMiranda & Soln 
(Il989h for the bright compact central region of the nebula. Our 
value for the Hubble-flow velocity scale factor was adjusted un- 
til the synthetic velocities from positions 4 and 5 matched the ob- 
served velocities; we retained a final value for the scale factor of 
9.5 km s^ Varcsec, and the resulting synthetic p-v arrays are shown 
in Fig. [2] (columns c and d). 

Our synthetic shell was tilted at various angles to the line of 
sight with the southern part of the shell towards the observer and 
the northern part away, and the results for /shell = 85°, 75°, 65° 
and 55° are displayed in Fig. |2] in rows 2, 3, 4 and 5, respec- 
tively. The model that best reproduces the observed p-v arrays is 
the /shell = 65° model (row 4). It becomes obvious from the syn- 
thetic p-v arrays in rows 2 and 3 that the PN shell is unlikely to 
have an inclination angle greater than 80° . 

As mentioned above, ^Kwok & Hsial ( |2007|) measured an in- 
clination of /shell = 52° from fitting el Upses to co-axial ring s ob- 
served in [Nil] HST images of Hb 12. iHvung & AUeJ ( 11996) cal- 
culate /shell = 40° from optical spectroscopy and photoionisation 
models. In t heir kinematical stud y of the nested double hourglass 
Hen 2-104, ICorradi et alj ( 1200 ll) found that the inner and outer 
hourglasses had identical inclination angles. We fitted inclined cir- 
cular rings to the equatorial features marked b in the infrared image 
of Kwok & Hsia ( 2007) and found an inclination of ~ 60° for the 
outer hourglass in Hb 12, which should be close, if not iden tical, to 
the inc lination of its bright inner hourglass. Finally. Hora & Latterl 
( Il996l) also found 60° from fitting an ellipse to the waist of the 
outer hourglass. 

Current theories of PN evolution suggest that a binary stellar 
system will shape the ejected PN shell such that the nebula shares a 
common inclination angle with the binary system. Hence an eclips- 
ing binary should produce a nebula with a symme try ax is at about 
90° (i.e. in the plane of the sky, see Mitchell et al] |2007i for exam- 
ple). All the measured values for the inclination of the hourglass 
shell of Hb 12 (65°, 60°, 52° and 40°), although not in strict agree- 
ment, together point to an inclination which is much lower than the 
80° limit which would be imposed if the central star of Hb 12 is in- 
deed an eclipsing binary system. This inconsistency in inclination 
angle could possibly be explained if the central syste m was a re- 
flection effect binary rather than a true eclipsing binary; iHsia et al.l 
( 12006) mention this as a possible interpretation of their lightcurve. 
One other way to maintain a common angle between the nebular 
shell and the progenitor system could be to invoke a triple system. 

As mentioned earlier, there is also some doubt concerning 
the eclipsing binary in Hb 12 (see iDe Marco et a l. 2008). In ad- 
dition to the arguments presented by De Marco et all the discrep- 
ancy between the binary and the nebular inclinations revealed by 
our study further suggests that the central stellar system is proba- 
bly not eclipsing. 
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Figure 2. Row 1 : optical data, (a) HST Ha image of the hourglass shell of Hb 12. (b) [Nil] A,6584 SPM image of Hb 12. (c) Position 4 [Nil] A,6584 p-v array, 
(d) Position 5 [Nil] A6584 p-v array. 

Rows 2-5 : Column (a): Raw shape particle model. Colvunn (b): raw model convolved with 2" seeing. Column (c): position 4 synthetic p-v array. Column 
(d): position 5 synthetic p-v array. Inclination angles are 85°, 75°, 65° and 55° for rows 2, 3, 4 and 5, respectively. The dashed lines in columns (a) and (b) 
mark the positions of sUts 4 and 5. 
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3.2 Evidence for bipolar ejection 

3.2.1 High-speed knots close to the central star 

In Fig. [3] we have displayed the Ha p-v arrays for slit positions 
2 to 6, which are all close to the central star of Hb 12, using a 
common velocity axis. For each frame, the slit is vertical, west is at 
the top and east at the bottom. Each frame measures 50" in size in 
the vertical direction. 

We can clearly identify a high-velocity ejection feature which 
is pointing in the direction away from the observer (positive veloc- 
ity) for slit positions north of the central star of the PN (positions 
5 and 6), and pointing towards the observer (negative velocity) for 
slits placed to its south (positions 2 to 4). This is the same north- 
south orientation as the PN shell. We also note that the velocity 
of the knots increase with distance from the central star. It is the 
first time that knots of this kind have been uncovered in Hb 12. 
I Welch et" al. ( 1999) detected [Fell] shock excited regions along the 
polar direction close to the core that should naturally form in the 
cooling wake of the bow-shock produced by a jet or bullet, thus 
predicting the presence of ejecta in or around Hb 12; our observa- 
tions confirm this for the first time. 

We fitted Gaussian line profiles to the emission features in or- 
der to obtain measurements for the velocities of the knots; an ex- 
ample of this is shown in Fig.|4] The measured knot velocities (cor- 
rected for systemic motion) and angular distances from the central 
star are listed in Table [T] By assuming that the inclination of ejec- 
tion direction is the same as that of the PN shell (i.e. i'ej = 65°) and 
that it runs down the polar axis of the hourglass, we computed the 
deprojected velocities at which the knots are moving away from the 
core, and their physical distances from the central star assuming a 
distance of 2.24 kpc. 

We have plotted in Fig.[5]the radial (systemic corrected) veloc- 
ities of the knots as a function of angular distance from the centre 
of the PN. A straight line through the origin 
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adequately fits the five inner data points (solid line), consistent with 
Hubble-type expansion. Hence, the observed ballistic ejecta most 
probably originate from a single eruptive event simi lar to the ejecta 
observed in NGC 6302 jMeaburn et alj2005Ll2008l) . 

We have also considered possible separate linear fits to the 
data points a bove and below the o rigin (dashed lines), similarly 
to the work of O'Connor et all |2000) on MyCn 1 8, although many 
less data points are available to us. Such separate fits would suggest 
that the ejection of the material would have taken place at some 
distance on either side of the central star It is interesting to note 
that this is the second hourglass PN for which this might be the 
case, even though no mechanism is currently known to produce 
such outflows. However, due to the uncertainties in the velocities 
and the low number of data points (for the points above the origin, 
we are 'fitting' a straight line to only two data points!), we conclude 
that the simple single straight line fit is a much better description 
of this central outflow. 

Finally, it is of course possible for the knots to have interacted 
with the circumstellar medium (CSM) which would most likely 
have caused them to decelerate, if the material is dense enough. 
The complex ejection of bursts of material at different velocities 
and/or different times into a CSM can only be properly treated us- 
ing hydrodynamical models, and we cannot predict here the be- 
haviour of the ejecta. However, we can show using a simple drag- 
force model (where the friction force acting on a moving clump is 





Velocity (km/s) 

Figure 3. Position-velocity arrays for slit positions 2 to 6. Velocities are he- 
liocentric. All spectra are Ha except for position 6 where the [Nil] A6584 
spectral line is shown, as we did not have an Ha spectrum at that slit po- 
sition in our dataset. West is up and east is down and each frame measures 
50" in height. 
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Table 1. Knot velocities, distances from centre and kinematical ages for the emission features observed in positions 1 to 7. The velocities have been corrected 
for systemic motion. Deprojected quantities and kinematical ages were calculated assuming i^j = 65° and D = 2.24 kpc. 



Slit 


Radial 


Deprojected 


Distance 


Deprojected 


Kinematical 


position 


velocity 


velocity 


from centre 


distance 


age 




Vr 




e 


d 


KA 




(kms-') 


(kms-') 


(arcsec) 


(xlO'^m) 


(years) 


1 


-125 ± 3 


-296 + 7 


72.3 + 0.3 


26.69 + 0.09 


2860 + 60 


2 


-75 ±4 


-178 + 10 


16.1 +0.3 


5.95 + 0.09 


1060 + 60 


3 


-48 ±4 


-114+ 10 


10.0 + 0.3 


3.70 + 0.09 


1030 + 90 


4 


-21 ±6 


-50 + 14 


3.7+0.3 


1.37 + 0.09 


870 + 230 


5 


+28 ±6 


+66 + 14 


4.2 + 0.3 


1.56 + 0.09 


750 + 100 


6 


+58 ±4 


+ 137 + 10 


11.5 +0.3 


4.25 + 0.09 


980 + 70 


7 


+ 116 + 5 


+275 + 13 


76.3 + 0.3 


28.16 + 0.09 


3250 + 90 
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Figure 4. Gaussian fits to the Ha emission line for position 3. 




Distance from centre (arcsec) 

Figure 5. Knot velocity as a function of observed angular distance from 
the central star Velocities are heliocentric and have been corrected for the 
systemic motion of the PN. The straight line fit to the five data points close 
to the origin follows « 6 . The dashed lines are separate fits to the data 
points above and below the origin. The two outer points correspond to the 
low-ionisation emission features at positions 1 and 7. 



proportional to its velocity) that if the CSM is uniform and all the 
ejected clumps have the same mass and different ejection veloci- 
ties, their position-velocity law would still go through the origin 
(Hubble-like, as in Fig. |5] solid line), thus rendering it indistin- 
guishable from a case with no friction using different initial veloci- 
ties. The gradient would only depend on the amount of friction (i.e. 
the density of the CSM) and the mass of the clumps. This scenario 



is somewhat idealistic and the probability for all these conditions 
to occur at once is rather low, which leads us to believe that the 
observed bullets would have only suffered a small amount of de- 
celeration. 

3.2.2 High-speed external knots 

Slit positions 1 and 7 were placed on knotty emission regions visi- 
ble in Fig.[T] and their [Nil] A6584 spectra (integrated over ~ 10" 
of slit length) are shown in Fig. [6] The radial velocities were mea- 
sured at Vr(l) = -125 ± 3 kms"' andV,(7) = 116±3kms-' for 
positions 1 and 7, respectively. These end caps are high-speed low- 
ionisation bipolar ejecta along the hourglass axis, mostly visible in 
the [Nil] A6584 emission line which implies coUisional ionisation 
by shocks, usually associated with regions where fast ejecta meet 
the interstellar medium (ISM) (see lLopez et al.|[l993l) . As for the 
inner knots, by assuming that the outer knots lie on the polar axis 
of the hourglass and that their ejection direction shares the same in- 
clination angle, we calculated their deprojected velocities and dis- 
tances; those are listed in Table[T] 

The compact nature of the knots suggest that they were also 
formed by bullet-like ejection. However, when plotted alongside 
the inner bullet velocities (see Fig. [5), the outer points lie very far 
from the Hubble-type expansion law. If they all originated from the 
same eruptive event (as seems to be the case for the knotty outflows 
in NGC 6302: iMeaburn et al.ll2005l |2008|) , the outer ejecta must 
have suffered severe deceleration (more than 50% of their initial 
velocity) from interaction with the ISM. This would also mean that 
the eruption ejected material on a very wide range of velocities 
K ~ 10 — 700 kms^' (700 km s^' would be the initial velocity of 
the outer knots if they followed the Hubble expansion law of the 
inner knots). The other possible explanations are that they belong 
to a different ejection event than the central ejecta and/or that they 
are at a different inclination angle. 

3.2.3 Kinematical ages 

The kinematical ages (KAs) for emission features I to 7 are listed 
in Table[T] The KAs for the end caps are approximately three times 
larger than for the inner emission features. It is thus likely that 
the inner and outer outflowing emitting gas belong to two sepa- 
rate eruptive events. The observation of separate eruptive events 
could be evidence that Hb 12 hosts a bipolar rotating episodi c jet 
or BRET ( iLopez, Meabum & PalniieJl 19931 : lL6Dez et aljl995l) . al- 
though it is not possible from our data to measure any rotation of 
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Figure 6. [Nil] A6584 spectra for slit positions 1 (top) and 7 (bottom) in- 
tegrated over 10" in slit length. The velocities have been corrected for sys- 
temic motion. 



the jet. One could also argue that the southern outer knot has a dif- 
ferent KA than the northern one indicating that themselves could 
have come from separate eruptive events. However, at such large 
distances from the central star, deceleration effects from interac- 
tions with the ISM are more likely to be significant, leading to dif- 
ferent values for KA for these knots even though they may have 
originated in the same event. 

Finally, we calculate from our synthetic model of the hour- 
glass nebular shell the KA of the nebula to be 1120 ±30 years, 
which is of the same order as the KAs of the inner knots. This sug- 
gests that the end caps were formed prior to the PN shell, as was the 
case in A63 (Mitchell et al. 2007), most probably during the end of 
the AGB phase or the proto-PN phase. There is much evidence in 
the literature for bipolarit y and ejection events in post-AGB stars 
(see lBalick & FraiiMl2002l : for example). 

This KA is muc h higher than the 300 years derived by 
iMiranda & Soi3 ( Il989l) . The difference probably comes from the 
fact that their KA only applies to the bright central region of Hb 12 
measuring 0.9" in diameter (which is not included in our model) 
compared to ours which applies to the entire bright hourglass. We 
suspect that these might in fact be two different components; the in- 
ner bright component may just be a bright thick stellar wind which 
is not kinematically part of the hourglass. Indeed, for a structure 
measuring 0.45" in radius, our hourglass Hubble law would give an 
expansion velocity of 4.3 km s^' compared to the 16 km s^' they 
used. Other differences might also come from the fact that the size 



of the component from which they derived their value was smaller 
than their seeing of 1.5", even though they corrected for instru- 
ment and seeing effects (in fact, they point out that 0.9" is an up- 
per limit and Sahai & Trauger 1998; measure it to be 0.6", which 
would bring the KA down to ~ 200 years for the same expansion 
velocity). 

After studying many young P Ne with point-symmetric mor- 
phologies, ISahai & Traugeil ( Il99 g) have suggested an interesting 
evolutionary scenario. A collimated outflow (or several of them) 
from the central star system is initially active during the late AGB 
stage and carves structures in the surrounding AGB wind. This 'sets 
the stage' for the development of an aspherical nebula where the 
hot fast wind (even if it is intrinsically spherical) runs into the im- 
printed AGB envelope. Such an evolution would be consistent with 
our calculated KAs. 



4 DISCUSSION AND CONCLUSIONS 

We have created a morpho-kinematical model of Hb 12 based on 
imaging and high-resolution longslit spectroscopy of the nebula. 
We found that the shell inclination at which the synthetic spec- 
tra best matched the observations was /shell = 65°, which is much 
lower than the inclination of a disputed eclipsing binary I'binaiy > 
80° derived from lightcurve studies. This result would go against 
current ideas of PN shaping from binary s ystems and th e com- 
mon binary-shell inclination as observed bv iMitchell et alj ( l2007h 
in A63, and supports the proposition that if Hb 12 does host a bi- 
nary central star, it is probably not eclipsing. 

Further analysis of the Ha longslit spectra revealed the pres- 
ence of a high-velocity bipolar outflows with Hubble-type veloci- 
ties close to the PN core oriented in the same north-south direction 
as the PN shell. High velocities were also measured in the longslit 
spectra of the end caps observed in the wide-field image of Hb 12, 
showing they were also associated with jet-like activity. It is the 
first time that such outflows have directly been observed in Hb 12. 
Even though they were oriented in the same direction, large differ- 
ences in velocities and KAs between the inner and outer knots led 
us to suspect that they originated from two separate eruptive events. 
We believe this may be evidence for a BRET in Hb 12. 

Binary systems are believed to be essential to the formation of 
BRETs, where the accretion disc around the primary feeds matter 
and energy in the high velocity outflows, and precession or oscilla- 
tions in the binary orbit cause the jets to precess likewise. This may 
suggest that Hb 1 2 does possess a non- eclipsing binary central star. 

As noted bv lKwok & HsialbOOTh , Hb 12 presents strong sim- 
ilarities to the hourglass nebulae Hen 2-104 and MyCn 18. The 
hourglass is first smooth close to the central star, and gets more ir- 
regular as we move away from the waist. The outer regions present 
a succession of rings, followed by some mottled structure where 
instabilities start to occur, with also some jet-like features emanat- 
ing from the rim. We have also presented in this paper evidence for 
knotty outflows along the polar axis, as is the case for Hen 2-104 
and MyCn 18. All three PNe present nested hourglass morpholo- 
gies and the 'eye' in the central regions of Hb 12 and MyCn 1 8 look 
very much alike. In the case of Hen 2-104, ICorradi et alj i200lh 
found that the inner and outer nebulae had similar KAs, suggest- 
ing simultaneous ejection. With a lack of spectroscopic data on the 
outer hourglass, we cannot verify this in the case of Hb 12. 

Hb 12 is the fourth object to have been discovered to have a 
nes ted hourglasses mor phology (Hen 2-104, MyCn 18 and R Aqr; 
see ISolf &UlriclJl985h . and in all three previous cases, the central 
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system has be e n sugg ested or confirmed to be a symbiotic binary. 
iKwok & Hsi3 ( 120071) proposed that the central system in Hb 12 
might also be a symbiotic binary (bearing in mind that there may 
be so me doubt on their derived orbital parameters : iDe Marco et alj 
I2OO8I) . The hourglass PNe could be uncovering a link between 
PNe, symbiotic stars and novae (see also the discussion bv lCorradil 
I2OO3I : on nebulae around symbiotic binaries). Moreover, this link 
could be strengthened by studies reporting novae outbursts oc- 
curing inside old bipolar PNe with somewhat narrow waists. A 
faint PN-like shell was detected around the remnant of GK Per 
iBode. O'Brien & Simps onll2004l) . and more recently around V458 
Vul ( Wesson et al., 2008) . 



ACKNOWLEDGEMENTS 

NMHV acknowledges his University of Manchester research stu- 
dentship during which most of the work presented here was per- 
formed, and his new ANR grant for the programme SiNeRGHy at 
the CEA/Saclay (ANR-06-CIS6-009-01). APR is supported by a 
STFC research studentship. JAL gratefully acknowledges financial 
support from DGAPA-UNAM grant INI 16908. We would like to 
thank the referee for sensible comments which have improved this 
paper. 



REFERENCES 

Balick B., Frank A., 2002, ARA&A, 40, 439 
Balick B., Preston H. L., Icke V., 1987, AJ, 94, 1641 
Balick B., Rugers M., Terzian Y., Chengalur J. N., 1993, ApJ, 411, 
778 

Bode M. E, O'Brien T. J., Simpson M., 2004, ApJL, 600, L63 
Bond H. E., 2000, ASPCS, 199, 115 

Cahn J. H., Kaler J. B., Stanghellini L., 1992, A&AS, 94, 399 
Ciardullo R., Bond H. E., Sipior M.S., Eullton L. K., Zhang C.-Y., 

Schaefer K. G., 1999, AJ, 118, 488 
Corradi R. L. M., 2003, ASP Conf. Ser., 303, 393 
Corradi R. L. M., Livio M., Balick B., Munari U., Schwarz H. E., 

2001, ApJ, 553,211 
De Marco O., 2006, lAUS, 234, 1 1 1 
De Marco O., 2009, ArXiv e-prints, 0902. 1 137 
De Marco O., Hillwig T. C., Smith A. J., 2008, AJ, 136, 323 
Dwarkadas V. V., Balick B., 1998, ApJ, 497, 267 
Erank A., 2006, lAUS, 234, 293 

Garcfa-Segura G., Langer N., Rozyczka M., Franco J., 1999, ApJ, 
517, 767 

Garci'a-Segura G., Lopez J. A., 2000, ApJ, 544, 336 

Gomez Y, Tafoya D., Anglada G., Miranda L. E, Torrelles J. M., 

Patel N. A., Hernandez R. E, 2009, ApJ, 695, 930 
Greaves J. S., 2002, A&AL, 392, LI 
Harrington J. P, Borkowski, K. J., 1994, BAAS, 26, 1469 
Hora J. L., Latter W. B., 1996, ApJ, 461, 288 
Hsia C. H., Ip W. H., Li J. Z., 2006, AJ, 131, 3040 
Hyung S., AUer L. H., 1996, MNIiAS, 278, 551 
Kahn E D., West K. A., 1985, MNRAS, 212, 837 
Kwok S., Hsia C. H., 2007, ApJ, 660, 341 

Livio M., Shankar A., Burkert A., Truran J. W., 1990, ApJ, 356, 
250 

Livio M., Soker N., 1988, ApJ, 329, 764 

Lloyd H. M., O'Brien T. J., Bode M. E, 1997, MNRAS, 284, 137 
Lopez J. A., Meabum J., Palmer J. W., 1993, ApJL, 415, L135 



Lopez J. A., Vazquez R., Rodriguez L. E, 1995, ApJL, 455, L63 

Meabum J., Lopez J. A., Gutierrez L., Quiroz E, Murillo J. M., 
Valdez J., Pedrayez M., 2003, Rev. Mex. de A. Ap., 39, 185 

Meabum J., Lopez J. A., Steffen W., Graham M. F. Holloway 
A. J., 2005, AJ, 130, 2303 

Meabum J., Lloyd M., Vaytet N. M. H., Lopez J. A., 2006, MN- 
RAS, 385, 269 

Miranda L. E, Solf J., 1989, A&A, 214, 353 

Miranda L. E, Gomez Y, Anglada G., Torrelles J. M., 2001, Na- 
ture, 414, 284 

Mitchell D. L., Pollacco D., O'Brien T. J., Bryce M., Lopez J. A., 
Meaburn J., Vaytet N. M. H., 2007, MNRAS, 374, 1404 

Moe M., De Marco O., 2006, ApJ, 650, 916 

Nordhaus J., Blackman E. G., Frank A., 2007, MNRAS, 376, 599 

O'Connor J. A., Redman M. P., Holloway A. J., Bryce M., Lopez 
J. A., Meaburn J., 2000, ApJ, 531, 336 

Porter J. M., O'Brien T. J., Bode M. E, 1998, MNRAS, 296, 943 

Sabin L., Zijlstra A. A., Greaves J. S., 2007, MNRAS, 376, 378 

Sahai R., Trauger J. T, 1998, AJ, 1 16, 1357 

Sahai R., et al., 1999, AJ, 118, 468 

Soker N., 1997, ApJ, 112, 487 

Soker N., 2001, ApJ, 558, 157 

Soker N., 2004, New Astronomy, 9, 399 

Soker N., Livio M., 1994, ApJ, 421, 219 

Soker N., Rappaport S., 2000, ApJ, 538, 241 

Solf J., Ulrich H., 1985, Ac&A, 148, 274 

Steffen W., Lopez J. A., 2006, Rev. Mex. de A. Ap., 26, 30 

Vlemmings W. H. T, Diamond P J., Imai H., 2006, Nature, 440, 
58 

Welch C. A., Frank A., Pipher J. L., Forrest W. J., Woodward 

C. E., 1999, ApJL, 522, L69 
Wesson R., et al., 2008, ApJL, 688, L21 



